Room temperature Peierls distortion in small diameter nanotubes 
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By means of ab initio simulations, we investigate the phonon band structure and electron-phonon 
coupling in small 4-A diameter nanotubes. We show that both the C(5,0) and C(3,3) tubes undergo 
above room temperature a Peierls transition mediated by an acoustical long-wavelength and an 
optical q — 2kF phonons respectively. In the armchair geometry, we verify that the electron-phonon 
coupling parameter A originates mainly from phonons ai q = 2kF and is strongly enhanced when 
the diameter decreases. These results question the origin of superconductivity in small diameter 
nanotubes. 

PACS numbers: 61.46.-l-w, 61.50.Ah, 63.20.Kr, 63.22.-fm, 68.35.Rh, 74.25. Kc, 74.70.Wz 



Superconductivity (SC) has been recently discovered 
in the 4-A diameter carbon nanotubes (CNTs) embed- 
ded in a zeolite matrix with a transition temperature 
Tsc= 15 K, much larger than that observed in the bun- 
dles of larger diameter tubes . This has stimulated a 
significant amount of work at the theoretical level to un- 
derstand the origin of the SC transition While 
both reduced screening and quantum/thermal fluctua- 
tions disfavor superconductivity in ID systems, early 
work on fuUerenes [jI and nanotubes suggested that 
the curvature significantly enhances the strength of the 
electron-phonon (e-ph) coupling. This is consistent with 
the softening of the Raman modes under tube diameter 
reduction 

In ID systems, the enhancement of the e-ph deforma- 
tion potential should also favor the occurrence of Peierls, 
or charge density wave (CDW), transitions driving the 
system into an insulating phase. This possibility has been 
explored in early works on the basis of model Hamilto- 
nians 01 ■ The strength of the e-ph coupling was ei- 
ther extrapolated from its value in graphite or taken as 
an adjustable parameter. The Peierls instability tran- 
sition temperatures Tcdw were estimated to be much 
smaller than room temperature: for instance, Huang et 
al. [13 find Tcdw 9.1 K for a C(5,5) tube. Further, 
long-wavelength optical Q or acoustical "opening 
the gap" phonon distortions have also been explored and 
the transition temperature was again predicted to be a 
few degrees Kelvin. 

The Peierls transition in 4-A diameter CNTs has 
not been addressed explicitly, except for a recent tight- 
binding exploration of the C(5,0) tube [l^ and an ab 
initio study of the C(3,3) one While the inclu- 

sion of electronic screening in the non-self-consistent ap- 
proach [i2 for the C(5,0) case yields a negHgible Tcdw, 
the ab initio simulations predict for the C(3,3) tube a 



Tcdw temperature of ~ 240 K. These results suggest 
that the C(5,0) tube is the only candidate for a SC tran- 
sition at 15 K. However, such a large discrepancy between 
the two results call for a unified treatment of both nan- 
otubes within a parameter-free approach. 

In the present work, we study by means of ab initio 
simulations the phonon band structure and e-ph cou- 
pling in 4-A diameter CNTs. We show that the C(5,0) 
and C(3,3) tubes both undergo a Peierls distortion with 
a critical temperature Tcdw larger than room temper- 
ature. While the distortion is associated with a long- 
wavelength acoustical mode in the C(5,0) case, the e-ph 
coupling is dominated by phonons a,t q = 2kp in the case 
of C(n,n) nanotubes with a coupling parameter A which 
is strongly enhanced with decreasing diameter. 

The electronic and vibrational properties are studied 
within a pseudopotential 0| planewave approach to den- 
sity functional theory (DFT). We adopt the local den- 
sity approximation (LDA) and a 50 Rydb cut-off to ex- 
pand the electronic eigenstates. The phonon modes are 
calculated at arbitrary g-vectors using the perturbative 
DFT approach [l^ (in what follows q- and fc-vectors la- 
bel phonon and electron momenta respectively). Special 
care is taken in sampling the electronic states around the 
Fermi level using an adaptative fc-grid technique. The en- 
ergy levels are populated using either a Gaussian broad- 
ening [3 or a Fermi-Dirac (FD) distribution within Mer- 
min's generalization of DFT to the canonical ensem- 
ble. The latter technique explicitly accounts for the effect 
of temperature. We use a periodic cell allowing for 10 A 
vacuum between the nanotubes. 

We start by studying the nanotubes band structure af- 
ter careful structural relaxation with a 25 meV FD distri- 
bution broadening (T'^300 K). The first important out- 
come of this calculation is the spontaneous zone-center 
deformation of the C(5,0) tube [l^l- From the Dioh sym- 
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FIG. 1: Symbolic representations of the out-of-plane acous- 
tical and optical modes at F for the phonon bands driving 
the Peierls instability in the (a) C(5,0) and (b) C(3,3) tubes, 
respectively. 




FIG. 2: Band structures of (a) the "distorted" D2h C(5,0) 
nanotube and (b) the C(3,3) tube. The zero of energy has 
been set to the top of the valence bands and at the Fermi 
level, respectively. 



metry, the structure relaxes to form an elliptic tube (D2h 
symmetry), as indicated by the arrows in Fig. ^a), with 
principal axes of 3.83 A and 4.13 A. 

The corresponding band structure is provided in 
Fig- EI a-) • As a result of the reduction of symmetry, the 
repulsion of the hybridized a-n bands at the Fermi 
level opens a ^ 0.2 eV band gap (LDA value). This 
implies that the C(5,0) nanotube undergoes a metallic- 
semiconducting transition with Tcdw larger than room 
temperature. Following Ref . |^ such a transformation can 
be referred to as a "pseudo-Peierls" transition since it 
opens the band gap through a (long-wavelength) struc- 
tural deformation. 

In the C(3,3) case, we do not observe any zone-center 
deformation leading the band structure away from the 
well-documented |20j zone-folding picture of two linear 
bands crossing at the Fermi level Ep- However, com- 
pared to the zone- folding analysis, the Fermi wavevec- 
tor fci? is no longer at 2/3 FX but at ~ 0.58 FX so that 
27: /2kp is no longer commensurate with the unperturbed 
lattice (or commensurate with a very large supercell). 
This is a well documented effect of the curvature. 

We further present in Fig. 01 a) the phonon band struc- 
ture of the C(3,3) tube 21]. The most salient feature is 




FIG. 3: Phonon band structures for the C(n,n) nanotubes 
with n=3, 4, and 5 from left to right obtained with a Gaus- 
sian broadening of 136 meV (10 mRy). The dotted line for 
the C(3,3) nanotube corresponds to the additional explicit 
calculation of the phonon modes around q = 2kF with a re- 
duced 68 meV (5 mRy) broadening. Important bands appear 
as thicker lines. The "undressed" ijj^ bands (see text) are 
indicated with thick dashed lines. 



the giant Kohn anomaly ^^25] at q — 2kp leading to a dra- 
matic softening of a few phonon modes |26| . In particu- 
lar, the optical band starting at 620 cm~^ at F becomes 
the lowest vibrational state at g = 2kp [thick/dotted 
low lying band in Fig.l^l^a)]. The associated atomic dis- 
placements, indicated in Fig.^b), correspond to an out- 
of-plane optical mode in agreement with the results of 
Ref.E This is consistent with the analysis provided by 
Mahan |2^: the coupling of electrons with out-of-plane 
optical modes should be enhanced away from zone-center 
and with decreasing diameter. The predominance of such 
a mode for e-ph coupling is a strong manifestation of the 
importance of curvature in small diameter tubes. 

The value of the associated frequency at q — 2k p 
strongly depends on the energy distribution broadening 
used in the phonon calculation, as clearly illustrated in 
Fig. Ola) by the difference between the thick and dot- 
ted hues. Using a 25 meV FD distribution (T-^300 K), 
we find a negative phonon mode u!2kp,u=i of -200 cm~^. 
This means that the C(3,3) tube undergoes a Peierls dis- 
tortion with Tcdw larger than room temperature '28]. 
Consequently, at the mean-field SC transition tempera- 
ture measured for the 4-A diameter CNTs (Tsc^l5 K), 
the C(5,0) and C(3,3) tubes should both be semiconduct- 
ing, seriously questioning the origin of superconductivity 
in small diameter tubes. 

Our results are in good agreement with those of Ref. ^3 
concerning the C(3,3) tube, but contrast significantly 
with those based on a model for the interaction Hamilto- 
nians and the electronic susceptibility We emphasize 
that our self-consistent treatment automatically includes 
the renormalization of both the c-ph and e-e interactions 
by electronic screening at the DFT level. 

To study the influence of tube diameter on the soften- 
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ing of the modes at g = 2kF, we represent in Fig. |3Ib,c) 
the phonon band structure of the C(n,n) tubes with n=4 
and 5 0| . Again, we observe that the same optical band 
(thick low lying band) is softened at 2kp , but clearly the 
effect is strongly reduced with increasing diameter. This 
confirms that Tcdw will quickly lower with increasing 
diameter. The large amount of sp^ character in ultra- 
small tubes invite to draw a comparison with the SC tran- 
sition in doped carbon clathrates and diamond [i^ . 

Surprisingly, the softening of another mode at higher 
energy (starting as an optical mode around 1450 cm~^, 
also indicated by a thick line) seems to be less sensitive 
to the diameter. This softening is the signature of a large 
e-ph coupling which might thus play an im por tant role 
in the transport properties of larger CNTs. [2J| 

To further quantify the strength of the e-ph interac- 
tion, we can calculate the e-ph coupling matrix elements 
gq„{knn') and the coupling constant A: 

X^Y,Xq^,^2N{Ef)Y,<\g,,\^> /nw,, (1) 



< \9qA > = 



L 



l^\9qu{knn)\ ^^^^^^ C^ ) 



BZ 



h 



SV 



CIV 



Due to conservation of energy and momentum, only vi- 
brational states arou nd o = T and ±2/eF contribute to A, 
as the nesting factor |23| n(q) quickly decays away from 
these points. At F, only four modes contribute to A due 
to symmetry considerations. They are the radial breath- 
ing mode (RBM), the optical out-of-plane (ZO) mode, 
and the in-plane optical longitudinal Ai(L) and trans- 
verse Ai(T) (G-band) modes (see Ref. 0) located respec- 
tively at 536 cm"-'^, 692 cm~^, 1354 cm~^ and 1431 cm~^ 
in the C(3,3) tube. 

In Table m we report the parameter \{i',q C {F}) de- 
fined as the g-sum of all the Xq,^, (see Eg. 1) for q in the 
neighborhood of F on a given z^-band [30j . This coupling 
strength varies from one band to another and is enhanced 
when the tube diameter decreases. However, the varia- 
tion with the diameter differs from one mode to another 
and it is difficult to extract a simple scaling law. 

We now discuss the e-ph coupling at g = 2fci?. 
The strong variation of the lowest frequency uj2kF,v=i 
with temperature raises the question of the meaning of 
^2kF,v=i calculated with the Qq^, and LOq^, values obtained 
from ah initio calculations. The difficulties have been 
clearly exposed in Ref. If one is interested in study- 
ing the SC transition, one faces the difficulty that at the 
experimental Tgc, the self-consistent g2kF,v and i02kF,v 
values are those of the CDW phase. As Tcdw is larger 
than any realistic SC transition temperature, we will 
not attempt to estimate Tsc with e-ph vertices and ph- 
propagator properly "undressed" from the CDW insta- 
bility jl2j . We limit ourselves to identifying the relevant 



TABLE I; Contributions (A^) to the coupling constant for all 
the g-vectors in the neighborhood of F (g C {F}) and 2fc_F 
(g C {2A;f}) for the C(n,n) tubes with n=3, 4, and 5. The 
corresponding frequencies {ijJv) at F and 2k f are indicated in 
cm~^. For the C(3,3) tube, the values underlined for g C 
{2^^} correspond to the "undressed" frequencies (see text). 
The global coupling constant (A) is indicated at the last line. 



C(3,3) 



C(4,4) 



C(5,5) 





UJi, 


A, 




A, 


UJi, 


A. 


gC{F} 
RBM 


536 


0.009 


416 


0.005 


338 


0.001 


ZO 


692 


0.009 


799 


0.004 


848 


0.001 


Ai(L) 


1354 


0.007 


1468 


0.009 


1538 


0.002 


Ai(T) 


1431 


0.012 


1505 


0.009 


1570 


0.002 


g c {2kF} 
















520 


0.074 


414 


0.022 


497 


0.010 




499 


0.026 


509 


0.010 


533 


0.004 




1320 


0.018 


1013 


0.010 


1154 


0.009 




1144 


0.001 


1180 


0.014 


1237 


0.011 






0.156 




0.083 




0.040 



modes and providing the values for X2kF,i'=i obtained 
above TcdWj which are those that matter for transport 
measurements in the normal state. 

While the g2kF parameters are rather stable with 
temperature above TcdWi the phonon frequencies still 
vary significantly. Following Ref. we provide an up- 
perbound for uj2kF,i' by removing the Kohn's anomaly. 
Namely, we assume that around q = 2k p, 



= K J' +Axln{q^ 2kF)/{q + 2kF) (4) 

and fit the bands by a cubic polynomial for [uj'^ ^)^ 
plus the logarithmic term that we subtract to obtain 
the bare phonon frequency. The important "undressed" 
bands (corresponding to w° j^) are represented by dashed 
lines in Fig. |3^a). The related frequencies and X2kF,v 
values are underlined in Tabled This treatment is ap- 
plied to both states presenting the largest Kohn anomaly, 
but its effect is clearly more pronounced for the low-lying 
V = 1 band. 



Using the u)'. 





2kF,v 



frequencies leads to providing a 



lower bound for the corresponding X2kF,i' e-ph param- 
eter. Even in that limit, we observe that the coupling 
with modes at g = 2kp is significantly larger than with 
long-wavelength phonons. As at zone-center, few modes 
with q — 2kF couple to the electrons and the strength 
of the coupling is significantly enhanced with decreasing 
diameter. In the C(3,3) case, we obtain in the undis- 
torted phase a density of states N{Ef) = 0.4 states/eV- 
cell-spin, yielding an e-ph potential — X/N{Ef) = 
440 meV. This is larger than that in the fuUerenes (60- 
70 meV), but it is not sufficient to lead to any signifi- 
cant Tsc value. The much larger density of states in the 
"Dioh" C(5,0) tube {N{Ef) = 1.8 states/eV-cell-spin) 
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may lead in principle to a large Tsc value. However, as 
shown above, the prevailing structure around the exper- 
imental Tsc temperature is the insulating one. 

In conclusion, we have shown that the "metallic" 
C(5,0) and C(3,3) nanotubes undergo a Peierls transition 
with a critical temperature larger than 300 K. Our re- 
sults are derived within Mermin's generalization of DFT 
to finite temperature systems and do not rely on models 
for the e-ph interactions nor the electronic susceptibil- 
ity. Additional physical ingredients (interaction with the 
zeolite network, possible defective structure, Luttinger 
liquid character, etc.) might reconcile these results with 
the experimentally observed SC transition at Tsc=15 K. 
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